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Tunneling conductance through two quantum dots, which are connected in series to left and right 
leads, is calculated by using the numerical renormalization group method. As the hopping between 
the dots increases from very small value, the following states continuously appear; (i) Kondo singlet 
state of each dot with its adjacent-site lead, (ii) singlet state between the local spins on the dots, 
and (iii) double occupancy in the bonding orbital of the two dots. The conductance shows peaks at 
the transition regions between these states. Especially, the peak at the boundary between (i) and 
(ii) has the unitarity limit value of 2e 2 /h because of coherent connection through the lead-dot-dot- 
lead. For the strongly correlated cases, the characteristic energy scale of the coherent peak shows 
anomalous decrease relating to the quantum critical transition known for the two-impurity Kondo 
effect. The two dots systems give the new realization of the two-impurity Kondo problem. 
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I. INTRODUCTION 



Dilute magnetic impurities in metal bring the single- 
impurity Kondo effect . The anti- ferromagnetic cou- 
pling between spins on impurities J, such as the RKKY 
interaction, would compete with the Kondo effect. To 
study such competition effect, the two impurities in metal 
have been studied extensively |^-[T^| . If the Kondo bind- 
ing energy is much larger than J (Xk S> J), each local 
spin on the magnetic impurity forms the Kondo singlet 
state with the conduction electrons. On the other hand 
for Tk <C J, the two local spins form the local spin sin- 
glet state. From the numerical renormalization group 
(NRG) calculation, Jones et al. had pointed out that the 
transition between the two states occurs as a quantum 
critical phenomenon p|. However the advanced investi- 
gation made clear that the critical transition is an artifice 
of the model neglecting the parity splitting terms, such 
as the d-d hopping term between the impurity atoms || . 
In this paper we will investigate the effect of the d-d 
hopping term in double-quantum-dot (DQD) systems in 
detail, and give the new realization to the two-impurity 
Kondo problem. 

It might be difficult to observe the two-impurity effect 
in metal systems as pointed out by previous studies, be- 
cause the alloy contains many types of impurity pairs, 
and because the coupling between impurities is fixed in 
each material. Recently, the Kondo effect is observed in 
the single quantum dot systems |l5|-[l^| . The experimen- 
tal data show good agreement with the results of numer- 
ical calculation based on the single-impurity Anderson 
model 



These works demonstrated that the quan- 
tum dot systems are suitable for sensitive experiment of 
the Kondo problem. On the DQD systems, each dot cor- 
responds to an impurity atom, and the coupling between 



the dots can be changed freely by applying the split gate 
voltage between the dots |2l[] . It would be expected that 
we can investigate the two-impurity effect systematically 
in the DQD systems. 

For the DQD systems, which the two dots are con- 
nected to the left lead and the right lead in series as 
'lead-dot-dot-lead', there are several theoretical works in- 
cluding the Kondo effect [2^[27[|. We have reported the 
large enhancement of the tunneling conductance through 
the two dots when the condition J£p ~ holds, by us- 
ing the NRG calculation (2f|. Here J£| = At 2 /U is the 
anti-ferromagnetic kinetic exchange coupling between the 
two dots, t is the hopping between the two dots, U is the 
Coulomb repulsion on the dot, and T K is the Kondo tem- 
perature at t = 0. We note that the anti- ferromagnetic 
coupling is the inevitable effect due to the kinetic pro- 
cess t and the Coulomb repulsion on the dot U. There 
are investigations with the slave boson mean field the- 
ory (SBMFT). Aono et al. had already studied the same 
model of us, however they could not find the relation 
Jlr ~ T K on the peak of the conductance pointed out 
by us because the SBMFT can not treat the kinetic ex- 
change process properly |24| ]. Georges et al. introduced 
the anti-ferromagnetic coupling J between the two dots 
by artifice in the model, and discussed the effect related 
to the critical transition on the conductance by using the 
SBMFT |M . However the introduction of the artificial 
,/ in the model and the calculation of the conductance 
within the SBMFT framework bring some questions as 
follows: (a) Does the effect of anti-ferromagnetic coupling 
pointed by Georges et al. actually appear in the DQD 
systems?, because the hopping itself breaks the quantum 
critical transition |^|5). If it appears, however, (b) how 
the conflicting effects of t, the kinetic exchange coupling 
that would cause the critical transition and the parity 
splitting that suppress the critical transition, compete? 
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And then, (c) how they appear in the conductance? Since 
the SBMFT could not treat the kinetic exchange process 
properly, this approximation for the two-impurity Kondo 
problem like the DQD systems seems to be unfavorable. 
The reliable calculation is necessary for such a sensitive 
problem. 

In this paper, we present the detailed investigation 
of the Kondo effect in the DQD systems. The numer- 
ical calculation is performed by using the NRG method. 
This numerical method is known to be a reliable one for 
the two-impurity Kondo problem |p|-p|,pp[p^] . We cal- 
culate the tunneling conductance through the two dots. 
We note that some preliminary results were presented 
at SCES98 fl2a], and the one of the central results was 
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presented at LT22 (27 1. 

We find that the following states continuously ap- 
pear when the hopping between the two dots increases 
from very small value; (i) Kondo singlet state (t -C U, 
J£ R « Tj), (ii) singlet state between local spins on the 
dots (t < U, J£ R > T£), and (iii) double occupancy m 
the bonding orbital of the two dots (i <; U). The conduc- 
tance shows peaks at the transition regions between these 
states. The 'main peak' at the boundary between (i) and 
(ii) with the condition J£ R ~ T^J has the unitarity limit 
value of 2e 2 /h because of coherent connection through 
the lead-dot-dot-lead. Especially for the strongly corre- 
lated cases, the width of the main peak becomes very 
narrow and the characteristic temperature of the peak 
is largely suppressed compared with the Kondo temper- 
ature of the single dot systems T£J. These anomalies of 
the main peak closely relate to the quantum critical phe- 
nomenon in the two-impurity Kondo problem. The quan- 
titative calculation in this paper gives the new realization 
for the two-impurity Kondo problem, and suggests the 
possibility of the systematic study of the anomalous two- 
impurity Kondo effect in the DQD systems. 

The formulation is presented in §[n| The numerical re- 
sults are presented in § III . The summary and discussion 



are given in § [V 



II. FORMULATION 

We investigate the following model Hamiltonian for the 
DQD systems that the two dots are connected to the left 
lead and the right lead in series; 

H = H\ + H d + fli_ d) (1) 
Hi = £kcl ka c Lka + Y £«4 gCT c Rg(T , (2) 
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+ (-tJ2 d LdR*+h.c.) 
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H\ gives the electrons in the left and right leads. gives 
that in the left and right dots. i?i_d gives the tunneling 
between the left lead and the left dot, and between the 
right lead and the right dot. The suffices L, R mean the 
left and the right, respectively. CLfco- is the annihilation 
operator of the electron in the left lead, rfur is that in 
the left dot. nd,Lcr = d^d^a is the number operator of 
the left dot. e k is the energy of the state k in the left 
lead. £d,L is the energy of the orbital in the left dot. The 
quantity t is the matrix element between the left and 
right dots, and we call it as the 'hopping' between the 
dots hereafter. C/l is the Coulomb interaction between 
the electrons in the left dot. Vl is the matrix element 
between the left dot and the left lead. 

Here we consider only the single orbital in each of the 
dots. This situation is justified in the case that the typ- 
ical energy splitting between the orbitals in the dot is 
larger than the typical broadening of the energy levels, 
fed A, and when the temperature is smaller than the 
typical Coulomb repulsion between the electrons in the 
dots, T < U @§§. (The Kondo effect is not impor- 
tant in the case of T <; U.) We consider only the on- 
site Coulomb interaction between the electrons. Further- 
more, we consider only the nearest neighboring tunnel- 
ing, between the dot and its adjacent-site lead, between 
the two dots. The energies £d,L, £d,R can be changed by 
applying the gate voltage on the dots. V^k (Vaq) can 
also be changed by applying the split gate voltage be- 
tween the left (right) dot and the left (right) lead, t can 
be changed by applying the split gate voltage between 
the left dot and the right dot. 

In this paper we consider only the symmetric case with 
respect to the exchange of the left and the right. This sit- 
uation is written with the following relations; £d = £d l = 
e d)R , U = U L = C/ R , and A = Al = Ar = tt|V|V' (A 
is the hybridization strength between the dots and the 
leads, V = V^k — Vnk, Pc is the density of states in the 
leads. Here we consider that there are no /c-dependence in 
the matrix element and the density of states.) The model 
can be mapped into the two-channel Anderson Hamilto- 
nian by the unitary transform for the operators of the 
dots and the operators of the leads |2f| . Furthermore we 
consider the situation that there is one electron in each 
dot by adjusting the gate voltage on the dot, (n<i,i,) = 1, 

(«d,R) = I- 

We solve the Hamiltonian by using the NRG method, 
and calculate the conductance from the current correla- 
tion function within the linear response theory |2^,^,^9| . 
(For detailed calculation of the conductance, see ap- 
pendix of Ref. @.) 

At zero temperature, the conductance can be re- 
written by using the effective parameters of the fixed 
point non-interacting Anderson Hamiltonian as follows; 

G = ^HAS c (0+)-AS o (0+)| 2 
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_ 2e 2 4(i cff /A off ) 2 

~ ~7T (1 + (t off /A eff ) 2 ) 2 ' (5) 

We have used the relation, Q p — z p /(uj — £p ff + «Ap ff ), 
z p = Af/A, (p = e,o) at T = 0. Here the suffix p de- 
note the even and odd parity orbitals in the two dots. We 
note that the even orbital is the bonding orbital, and the 
odd orbital is the anti-bonding orbital. We now consider 
the case of (n e ) + (n Q ) = 2, then t cS = -ef = ef > 0, 
A eff = Af = Af. Here t cS is the effective hopping 
between the dots, and A eff is the effective hybridization 
strength between the leads and the dots. At T — 0, we 
calculate the effective parameters from the analysis of the 
flow chart of the renormalized energy level structure in 
the NRG calculation, and then calculate the conductance 
from eq. 
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FIG. 1. Conductance G, occupation number in the odd 
orbital of the two dots (n ), and phase shift of the odd chan- 
nel So, as a function of the hopping t at T = 0. We note the 
relations between the even and odd orbitals, (n c ) = 2 — (n ), 
S c = ty - 5 . A/tyU = 1.5 x 10" 2 . 



III. NUMERICAL RESULTS 



In numerical calculation we choose the half of the band 
width as an energy unit. The Coulomb repulsion is fixed 
at U = 0.1 through this paper. We calculate the conduc- 
tance as a function of the hopping t for various hybridiza- 
tion strength A. (As noted in previously, t and A can be 
changed by applying the split gate voltage between the 
dots, between the dots and the leads, respectively.) The 
gate voltage on the dots is fixed at £d = —U/2, then the 
DQD is in the half-filled case, i. e. each dot contains one 
electron. 
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FIG. 2. Density of states on the even orbital of the 
two dots, p e {io). We note the relation p (u) = p c (—La), 
where p (u) is the density of states on the odd orbital. 
A/tyU = 1.5 x ltr 2 . 



A. Conductance in the strongly correlated case 



First we present the conductance in the strongly cor- 
related case with the hybridization strength satisfying 
A/tt = 1.5xl(T 3 , (A/nU = 1.5xlCT 2 , i.e. u = U/ttA ~ 
6.8). 

We show the conductance at T — as a function of the 
hopping t in Fig. |l|. (The occupation number and the 
phase shift are also shown in Fig. |l].) There are two peaks 
in the conductance, the large peak near t ~ 5 x 10~ 4 , and 
the small peak near i~2x 10 -2 . (Hereafter we call the 
large peak as the 'main peak'.) Why these peaks appear? 
In the later paragraph we will analyze various quantities 
for the parameter cases showing the peaks. 



The density of states on the even orbital of the two 
dots, p e (uj), for several t cases is shown in Fig. ||. (The 
relation p c {uj) — p {—uj) holds for (n c ) + (n Q ) = 2, where 
/0 o (w) is the density of states on the odd orbital of the two 
dots.) At t — 0, there is the Kondo peak on the Fermi 
energy. (Fermi energy corresponds to u) = 0.) Naturally, 
this Kondo peak is caused by the Kondo singlet states be- 
tween the left lead and the left dot, and between the right 
lead and the right dot. As t increases to t = 5 x 10~ 4 , 
the conductance has the main peak and the strength of 
p e (u> ~ 0) becomes half of that at t = 0. In this region 
we can consider that the Kondo effect with the spins on 
the orbitals extending the two dots, the even and the 
odd orbitals, occurs. As t still increases, the conduc- 
tance decreases rapidly, and the strength of p c (ui ~ 0) is 
largely suppressed as shown att=1.0xl0~ 3 . This sup- 
pression means the disappearance of the Kondo coupling 
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between the leads and the dots, t becomes still large, 
the conductance has the small peak at t ~ 2 x 10~ 2 . 
At t = 2.5 x 10~ 2 , the density of states on the even and 
odd orbitals have peaks at =Fw ~ 10 -1 , respectively, from 
Fig. ^| At the same time the occupation numbers begin 
to change as shown in Fig. [y. ((n c ) ~ 1.5, (n Q ) ~ 0.5 at 
t~ 1.5xl0- 2 .) 

Here we note the following two points: First, the condi- 
tion J£g ~ T° holds at t ~ 5 x 10~ 4 , where J£^ = 4t 2 /U. 
(T$ = 3.78 x 10~ 6 is the Kondo temperature at t = 0, 
with the expression T£ = ^UA/2 cxp[-7rE//8A + A/2C7] 
§, then J£r/T5_~ 2.65 at i = 5.0 x 10~ 4 .) Second, as 
seen from Fig. |lj, the occupation numbers of the even 
and the odd orbitals for t <, 1.0 x 10~ 3 arc almost same 
with each other, (n c ) ~ (n ) ~ 1. For t ^ 1.0 x 10 _1 , 
the two electrons occupy the even orbital. The border 
between them is at t ~ U/4(= 2.5 x 10~ 2 ). 

Above analysis implies the following scenario. In the 
case of t <C U/4(= 2.5 x 10~ 2 ), the hopping t causes 
the anti- ferromagnetic kinetic exchange coupling, J£^. 
For smaller hopping case with J££ < T°, there are the 
Kondo singlet states between the left lead and the left 
dot, and between the right lead and the right dot each 
other. As t increases and then J£^ 3> T^, the two lo- 
cal spins on each dot form the local singlet state. At 
the transition region between two states we have a main 
peak with the unitarity limit value of 2e 2 /h in the con- 
ductance. This will indicate that the leads and the dots 
are coherently connected by the even and the odd or- 
bital states. When t becomes still large and the condition 
t <; U/4 holds, the local spins do not appear, instead, the 
two electrons occupy the even orbital. The small peak of 
the conductance reflects the transition of the electronic 
states in the DQD. 
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FIG. 3. Effective parameters t , A of the fixed point 
non-interacting Anderson Hamiltonian, given by the analysis 
of the flow chart of the renormalized energy level structure in 
the NRG calculation. A/tvU = 1.5 x 10~ 2 . 

Here we show the effective parameters t eS and A eS as 
a function of t in Fig. |[ (We note that the effective 
parameters have been used already for the calculation of 



the conductance shown in Fig. |].) In t < 10 4 case, the 
effective parameters behave as t cS ~ tlQ/A, A cS ~ Tg. 
Then the conductance coincides with the non-interacting 
one when one substitutes the effective parameters into eq. 
(||). As t increases with the condition near J£^ ~ T£, 
A eff once decreases, and it has local minimum, and then 
it increases. At the same time the slope of t once 
decreases and then increases. When t cS and A eff co- 
incide with each other, the conductance has a peak at 
(~5x 10~ 4 , i. e. J£p ~ T£. As t increases slightly be- 
yond this point, the conductance sharply decreases be- 
cause A eff decreases to the minimum even though t eS 
increases. Here we stress that the relation J£p ~ T£ 
holds when t cff ~ A cff , and at the same time A eff be- 
comes very small in the transition region. When t in- 
creases further, the ratio t /A increases gradually in 
the region t < £7/4(= 2.5 x 10~ 2 ). At t - £7/4, the ratio 
i cff /A cff begins to decreases and then increases. There- 
fore the conductance shows a broad peak near the region 
t ~ U/4. For t <; J7/4 case, the effective parameters 
behave t cS ~ t, A cS ~ A. We note that the conduc- 
tance has the expression of the non-interacting one itself 
in t ^ U/4 region. 

Finally we compare between the phase shift and the 
occupation number shown in Fig. [l]. The phase shift of 
the odd orbital 6 a rapidly changes from 7r/2 to near 
f~5x 10 -4 , even though (n ) still remains at (n Q ) ~ 1. 
Friedel's sum rule in each channel does not hold, as al- 
ready pointed out previously 0]. It seems that this be- 
havior is enhanced when the anti-ferromagnetic coupling 
between the two sites competes with the Kondo effect. 



B. From Weakly to Strongly correlated cases 



In this subsection we present the numerical results of 
the conductance for various A/irU cases within 1.5 x 
10- 2 < A/ttC7 < 6.0 x 10~ 2 . (1.7 <S u <, 6.8. The hy- 
bridization strength is changed in 1.5 x 10~ 3 < A/7T < 
6.0 x 10~ 3 , and the Coulomb repulsion is fixed at U = 
0.1.) We confirm the scenario shown in the previous sub- 
section that J£^ ~ T£! holds at the main peak of the 
conductance with t ~ A e . We also demonstrate how 
the kinetic exchange process appear in the conductance 
for arbitrary A/ttU cases. 

The calculated conductance is shown in Fig. [|. The 
horizontal axis is the hopping normalized by the hy- 
bridization strength, t/A. From inset of the Fig. 
the conductance almost overlaps on the non-interacting 
curve in the region t <C A and t ^> A. Exactly, these 
regions should be classified J£§ <§; 7^ and t 3> U/4, re- 
spectively, from the analysis in the previous subsection. 
The conductance is very small in these regions, however 
this uniform properties should be useful to arrange the 
experimental data under uncertain U/A cases. 
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FIG. 4. Conductance as a function of t at zero tem- 
perature, from the weakly to strongly correlated cases in 
1.5 x 1(T 2 < A/tyU < 6.0 x 10~ 2 . The broken line shows 
the conductance for the non-interacting (U=0) case. 

All curves have a main peak with strength 2e 2 /h. For 
weakly correlated cases of A/ttII ^ 4 x 10~ 2 , the con- 
ductance almost coincides with the non-interacting one 
through all t region. As U/A increases, the main peak 
shifts to the small t/A side, and the peak width becomes 
narrower. 

We already found the relation J£^ ~ T® at the main 



case in 



III A 



peak position for A/ttII = 1.5 x 10 
Here we show the ratio J£r/^k a ^ * ne mam peak po- 
sition for various A/ttU cases in Table j. We can see 
T£ commonly |25|. (The rela- 



the relation J£p 
tion at the main peak, 



Tcff 

J LR. 



Tjf, would be general- 



ized to Eb ~ ?k ^0 the weakly correlated cases, where 
E B = y/(2t) 2 + {U/2) 2 - U/2 is the singlet binding en- 
ergy between the two dots.) 

From the analysis in the previous and present subsec- 
tion we can conclude the following effect of the hopping 
term. For the small t case with J£| < T°, "(i) The 
Kondo singlet state is formed on the left (right) dot with 



its adjacent-site lead". On the other hand for the large 
t case with J£^ 3> T-£, "(ii) the local spins on each of 
the dots couple as the singlet state" . In the intermediate 
region, the Kondo effect of the local spins on the orbitals 
extending on the two dots (i.e. even and odd orbitals) oc- 
curs. The main peak of the conductance appears around 
the boundary between (i) and (ii) reflecting the coherent 
connection of the leads and the dots. As U /A increases, 
the Kondo temperature T£ exponentially decreases, the 
condition J£^ ~ holds at the smaller t/A, then the 
main peak shifts to the smaller t/A side. At the same 
time, the width of the peak becomes extremely narrow 
compared with the decreasing of T£. This fact has been 
already shown as the steep minimum of A cff in Fig. §. 
We note that this narrowing closely relates to the quan- 
tum critical transition between the Kondo singlet state 
and the local singlet state in the two-impurity Kondo 
model Q. The shifting and narrowing behaviors shown 
here are also pointed out with the SBMFT with artifi- 
cial addition of the anti-ferromagnetic coupling between 
dots to the model However, the SBMFT calculation 
should be checked by the method treating the kinetic ex- 
change term properly. As noted in the introduction, the 
hopping term causes two conflicting effects on the crit- 
ical transition of the two- impurity systems. One is the 
kinetic exchange coupling J£^, which causes the "criti- 
cal" transition through the competition with the Kondo 
effect. Another is the parity splitting, which suppresses 
the "critical" transition. The calculation in this section is 
the first reliable quantitative results of the two-impurity 
Kondo problem in the DQD systems. 

There is also another small peak (or shoulder) struc- 
ture for the strongly correlated cases of A/ nil ^2x 10~ 2 
(u > 5) at larger t side of the main peak. In the previ- 
ous subsection, we found that the small peak appears 
around the boundary between (ii) and (iii). However for 
the weakly correlated cases, the small peak could not 
be recognized because the condition of the border (i)-(ii) 
and (ii)-(iii) could not be distinguished clearly. 



A/tvU 


1.5 x 10" 2 


2 x 10" 2 


3 x 1(T 2 


4 x 1CT 2 


6 x 10 -2 


tCH /rpU 

Eb/T£ 


2.66 
2.66 


2.34 
2.34 


2.15 
2.14 


2.09 
2.04 


2.23 
2.05 



TABLE I. Ratios J£r/t£ and Es/Tk at the main peak position of the conductance. 



C. Temperature dependence of the conductance 



h to 



(6) 



In this subsection we present the conductance in fi- 
nite temperature. We calculate the conductance at finite 
temperatures by using the following formula 129]; 



Here P (uS) is the 'current spectrum' for the current op- 
erator J = TVl — N-r written as follows, 
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4 Z 



x\(n\J\m)\< 
x<5 (w - (E n 



-0E„ 



E rn )) 



-0E n 



(J) 



where Nt, is the time differentiation of the electron num- 
ber in the left lead, Z = J2 n e~^ En is the partition func- 
tion of the system, and (3 is the inverse of the temperature 
09 = 1/T). 

First we show the conductance at various temperatures 
for A/ttII = 1.5 x 10~ 2 case in Fig. [| As the temper- 
ature increases from T ~ 10 -8 , the height of the main 
peak gradually decreases. At the same time the peak po- 



sition shifts to the larger t. We note that T 
much lower than T°. (T° = 3.78 x 1(T 6 .) 
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FIG. 5. Temperature dependence of the conductance. 
Main figure is the temperature dependence near the main 
peak at t ~ 5 x 10~ 4 . Inset figure at the upper right is the 
temperature dependence near the small peak at t ~ 2.5 x 10 -2 . 
Inset figure at the lower right is the conductance in all over t. 
A/tyU = 1.5 x 1(T 2 . 

To discuss the characteristic behaviors of the conduc- 
tance in finite temperature, we show the density of states 
Pc(uj), and the current spectrum (divided by uj) P [u))/w 
at t = 5.0 x 10~ 4 in Fig. ^|. As the temperature in- 
creases to T = 6.44 x 10~ 8 , P (u>)/u> at ui ~ become 
60% of T = limit. At the same time, p e (oj) shows a 
small change around ui ~ 10~ 7 . This means that the 
effect of the temperature on the conductance is rather 
drastic. Here we show two sorts of the magnetic exci- 
tation spectra x m (w) and f|' wnere Xm( w ) i s the 
imaginary part of the dynamical susceptibility of the uni- 
form magnetic moment of local spins, (Sl,z + Sr^Vv^, 
and x a (w) is that of the anti-ferromagnetic moment, 
(5l, 2 — Sr iZ )/V%, respectively. We show the two mag- 
netic excitation spectra at t — and t — 5.0 x 10~ 4 in 
Fig. [j| At t = 0, the two spectra agree with each other. 
However at t = 5.0 x 10~ 4 , Xa( w ) the structure in 



lower energy region than x m (w). It seems that P (u>) at 
the main peak of the conductance is dominated by the 
fluctuation given by x a ( w ) from Fig. |^ and Fig. ^. 
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FIG. 6. Density of states p B (u>) and current spectrum 
P (w)/w at t — 5.0 x 10 -4 for finite temperatures. (The 
spectrum, which has only the positive region u > 0, is the 
current spectrum.) 
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FIG. 7. Two sorts of the magnetic excitation spectra 
X m (w) and Xa ( w ) at T = 0. We note that the two spectra 
agree with each other at t — 0. 

We determine the two characteristic energies from 
X m (w) and Xa( w ) as the following ways. One is deter- 
mined from the peak position of Xm( w )> we can ^ 
pfj| , p8 129 1 . Another one, we call it Taf, is determined 
as T AF /T AFj0 = X /X, where X = lim^ \' a (w)/cj [§. 
(And here we have Taf.o = 2m,o-) The suffix '0' indicate 
H = 0'. The quantity Tm,o almost coincides with TS. 
The ratio Tm,o/^k f° r some cases are shown in Ref. [ po[ . 

The calculated two characteristic temperatures, Tm 
and Taf, are shown in Fig. ||. They take almost same 
values in t < 10~ 4 . Tm monotonically increases as t 
increasing. On the other hand Taf becomes smaller 
near J£ R ~ T^. It has minimum of Taf — 3 x 10~ 7 
at t ~ 5 x 10~ 4 . We note that the reduction of Taf 
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near T° 



J had been already pointed out Iq). As t 
still increases, Taf rapidly increases. From same analy- 
sis for the other A/irU cases, we confirm that the min- 
imum of Taf appears for the strongly correlated cases 
of A/ttU < 2 x 1CT 2 . We show Taf at the main peak 
position in Table. ||. 




FIG. 8. Two characteristic energies Tm and Taf- 

From the comparison with effective parameters in Fig. 
||, the larger of the effective parameters, max(£ eff , A eff ), 
and the smaller of the characteristic temperature, 
min(TAF)TM), almost coincide with each other in all t, 



max(t cB , A cff ) - min(TAF,T M ). 

Here we again see the temperature dependence of the 
conductance shown in Fig. |^ with the characteristic tem- 
perature shown in Fig. |§|. We can see that Taf charac- 
terizes the main peak of the conductance in finite tem- 
perature. The peak decreases as the temperature in- 
creases near T ~ 1 x 10 _8 (~ 0.1T AF (i = 5 x 1CT 4 )) 
in Fig. ^ As the temperature increases and reaches to 
T~ lxl(T 6 (~ 10T AF (t = 5xl(T 4 )), the strength of the 
main peak becomes almost zero. Next we see the tem- 
perature dependence of the small peak. The small peak 
near t ~ U/4 = 2.5 x 10~ 2 increases as the temperature 
increases to about T - 1(T 3 (~ 0.1T M (t = 2.5 x 1CT 2 )). 
It seems that the characteristic temperature of the con- 
ductance near the small peak is Tm- From above it seems 
that the characteristic temperature of the conductance is 
min(TAF, Tm) in all t region. 

Here we show the conductance from the weakly to 
strongly correlated cases at fixed temperatures. We show 
the conductance at T = 1.6 X 1CT 5 and T = 1.4 x 1CT 4 in 
Fig. ^. The main peak for the strongly correlated cases 
is sensitive to the temperature. Then the main peak of 
the conductance will shift to smaller t/A side with in- 
creasing peak height when the temperature decreases as 
seen from Fig. ^. This behavior will be observed as the 
split gate voltage is varied. We note that T — 1.4 x 10~ 4 
corresponds to 16mK, and T = 1.6 x 10~ 5 corresponds 
to 1.9mK, for U = l.OmeV systems. 



A/nU 


i.5 x lcr 2 


2 x 1CT 2 


3 x 1(T 2 


4 x 1(T 2 


6 x 10~ 2 


Taf 


3.28 x 1(T ' 


9.77 x 10"° 


3.97 x 1(T 4 


1.50 x 10* 


3.01 x 10* 


Tm 


7.66 x 1(T 6 


5.64 x 1(T 5 


7.73 x 1(T 4 


3.47 x 10 -3 


1.20 x 10 -2 


( Tm,o 


2.55 x 1(T 6 


2.47 x 10~ 5 


2.31 x 1(T 4 


8.02 x 10~ 4 


3.22 x 10~ 3 ) 



TABLE II. Characteristic energies at the main peak position. 
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Finally, we note the accuracy of the conductance cal- 
culated from eqs. (§)-(0) by using the NRG method. 
It is not so accurate at very high temperatures for the 
small peak. In the case of t = 0, two dots completely 
decouple, then the conductance should be zero. However 
as shows in Fig. |To|, the calculated conductance has fi- 
nite value in 5 x 10~ 3 <, T <, 1 and it has maximum at 
T ~ 0.05(= U /2). Thus the result at very high tempera- 
tures has ambiguities. This improper finite conductance 
would be caused by estimation of eq. (^J) at to ~ T in- 
stead of lu ~ * 0. The finite value in the current spectrum 
at T ~ 0.05 would reflect the largeness of the dynamical 
charge fluctuation in the dots. 



0.10 



0.08 



0.06 



O 0.04 



0.02 




'l0~ 10 10~ 9 10" 8 10" 7 10" 6 10" 5 10" 4 10" 3 10" 2 10" 1 1 
T 

FIG. 10. Numerical results of the conductance for t = 



peak became steeply narrow. The characteristic temper- 
ature of the main peak was strongly reduced compared 
with the Kondo temperature of the single dot systems 
These anomalous behaviors of the main peak related 
to the quantum critical transition of the two-impurity 
Kondo problem studied in previously. Though the hop- 
ping term had conflicting effects on the critical transi- 
tion of the two-impurity Kondo systems, generation of it 
through the kinetic exchange coupling and suppression 
of it due to the parity splitting, we found that we see the 
sign of the anomaly in the tunneling conductance. 

The quantitative calculation shown in this paper gave 
the new realization for the two-impurity Kondo problem. 
This investigation suggested the importance of the sys- 
tematic study of the DQD systems for the two-impurity 
Kondo problem. 
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IV. SUMMARY AND DISCUSSION 



We calculated the tunneling conductance through the 
two quantum dots that connected to the left lead and 
the right lead in series. We investigated the effect of the 
kinetic exchange coupling between the dots, and also the 
competition with the Kondo effect. 

As the hopping between the two dots increased, (i) 
Kondo singlet state, (ii) local spin singlet state, and (iii) 
molecular orbital like state with double occupancy in 
even state, continuously appeared. For t -C U cases, 
the Kondo binding between the left (right) lead and the 
left (right) dot, T^, and the anti- ferromagnetic kinetic 
exchange coupling between the two dots, J£^, competed. 
The boundary between (i) and (ii) was characterized as 
Jlr ~ TS, and the tunneling conductance showed a peak. 
This peak had the unitarity limit value of 2e 2 /h reflecting 
the coherent connection through the lead-dot-dot-lead. 
At t ~ f//4 of the boundary between (ii) and (iii) , we 
had a small peak. 

The system showed the strongly correlated behaviors 
for A/ nil & 2x 10" 2 (u = U/nA £ 5) cases. The borders 
of (i)-(ii) (J£| - T°) and (ii)-(iii) {t ~ U/A) were clearly 
distinguished, then there were the two peak structures 
in the conductance. Furthermore the width of the main 
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